The comparison between observations of galaxy clusters thermo-dynamical properties and theoretical predictions suggests that non-gravitational heating needs to be added into the models. We implement an internally self-consistent heating scheme into GADGET-2 for the third (and fourth) run of the Millennium gas project (Pearce et al. in preparation), a set of four hydrodynamical cosmological simulations with N=2(5x10^8) particles and with the same volume (L=500 h-1 Mpc) and structures as the the N-body Millennium Simulation (Springel et al. 2005). Our aim is to reproduce the observed thermo-dynamical properties of galaxy clusters.